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' Abstract 

Results of perturbation theory in quantum field theory generally de- 
' pend on the renormalization scheme that is in use. In particular, they 

^ I depend on the scale. We try to make perturbation theory scheme invari- 

■ ant by re-expanding with respect to a scheme invariant quantity. Further- 

more, we investigate whether the potentially large logarithms in such an 
, expansion cause inaccuracy and how this can be improved. 

^ ! 

^ ; 1 Introduction 

' The occurrence of divergencies in perturbative quantum field theory has made 

it clear that although measurable quantities should be finite, this need not 
^ i' be the case for the parameters of the theory. To handle this, divergencies are 

O . regularized. To obtain an order-by-order finite perturbation expansion, it should 

be specified what to do with divergencies as they occur at each loop level. 
This introduces an arbitrary constant for each loop level. Although the full 
perturbation series should not depend on these constants, the truncated one 
does. This is a problem if one is trying to approximate a physical quantity that, 
' by the definition of "physical quantity" , should not depend on these unphysical 

parameters. 

It has been pointed out in |1(J| that scheme invariant quantities can be con- 
structed out of the scheme-dependent ones. Therefore, the natural thing to do 
would be to try to rewrite the perturbation series into a series expansion with 
respect to the quantity Xi , the invariant that can be calculated at the one loop 
level. Actually, that would be a series expansion in l/Xi, since 1/A"i is the 
small quantity if we are in the perturbative regime. 

An expansion in 1/Xi will, however, contain logarithms in the expansion 
coefficients. At high energies, these logarithms will dominate all other con- 
tributions to the expansion coefficients. We will resum these logarithms and 
investigate whether this gives more reliable results. 

*C. Dams@science.ru.nl 
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X 



1 



Several attempts to remedy this problem of renormalization scheme depen- 
dence have been proposed. They all suffer from a great deal of arbitrariness 
and/or the mathematics involved is more complicated than that of simply ma- 
nipulating power series, as we will be doing. 

To obtain our results, we used the simple case of a one-parameter theory. 
We did not yet consider more complicated cases, but will do this later. When 
doing explicit calculations, we use QCD with five massless quarks, as a concrete 
example of a one-parameter theory. 



2 Renormalization Scheme Invariants 

In this section we outline how renormalization scheme-independent quantities 
can be combined into scheme invariant ones. This is also explained in | 10) and 
[Hj. For self-containedness of this paper and also to make clear what our con- 
ventions are, we will repeat this in this section. The coupling constant will 
be denoted by a. The key idea is that the consistency of perturbation theory 
requires that a result up to terms of order a" should only differ from the exact 
answer up to terms of order a""*"^. If we consider the case n = 1, we see that the 
running of the coupling makes the use of a as an expansion parameter scheme 
dependent. This indeed is an effect that starts at order a^. I.e., the lowest order 
term in the beta function is of order a^. 

The renormalization scheme can be specified by giving the renormalization 
scale s and the scheme-dependent beta function coefhcients f32, Pa, Pa, ■ ■ ■ ■ The 
coupling constant depends on the scale via 

^=P{a)^Poa^+Pia^+p2a^ + ---, (1) 

where ds = dfi/ fi with fi the mass scale used in dimensional regularization. This 
differential equation can be integrated to give 



'/3oa /?g Pia V/3(«') Poia'^ + Piia'f ' ' ^ ' 



This solution implies the choice of a boundary value. Our choice is similar to 
the one in |TQ]. Taking the set of variables that consists of s and the scheme- 
dependent beta function coefficients as independent, we can derive that 
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As an example, we imagine that we have calculated a physical quantity R 
up to fourth order. From this example, we hope, it will be clear how this can be 
generalized to arbitrary order. Having calculated R up to fourth order means 
that we have 

R r^a + ria^ + r2a^ + r^a^ . (4) 



Consistency of perturbation theory requires that 



To be more concrete, independence of s implies that 

(ro + 2ria + Sraa^ + 4rsa^) — + —^a^ + —^a^ + —^a^ ~ a^. 6 

^ ' OS OS OS OS 

From this, equations for second, third and fourth order can be extracted. Sub- 
stituting the beta function for da/ds, we find 

roPo + ^ = 0; 
2/?ori+ro/?i + ^ = 0; (7) 

OS 

3(^or2 + 2/3iri + ro02 + ^ = 0. 

OS 

To obtain the equations that follow from independence of (32, we expand da/ 8^2 
up to fourth order. We have 

Using this and then demanding that the third and fourth order of dR/d^2 are 
zero, we get the equations 

(3o dP2 ' 



2!l + ^ = o. 

/?o dp2 

Finally, independence of gives the equation 

^0 , drs 



(9) 



Integrating the equations for n, r2, and rs, we find 
ri=Xi- (3oros; 



r-2 = ^2 



ro /?o /3o ' (11) 



^3 = ^3 H H 



ro 2/?oro 2/3o /3o r^' 

where the quantities arc the renormalization scheme invariants. They arise 
because the values of the rj are obtained from differential equations and hence 
need constants of integration. The values of the Xi can bo calculated from the 
values of and j3i once these have been calculated. It should be made sure 
that the rj and Pi have been calculated in the same renormalization scheme, of 
course, and then it will turn out that the values of Xi no longer depend on the 
scheme that was used to obtain rj and Pi. 

The acute reader will have noticed that we actually expressed r2 and rs 
into ri instead of in s. Hence, we actually label our renormalization scheme by 
I"!, P2, Pz, ■ ■ ■ rather than by the set of parameters that we mentioned earlier. 
The reason that we do this is that if we are going to re-expand into 1/Xi we 
would rather have ri than s in the expansion coefficients because in a suitable 
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renormalization scheme ri will be a number of order unity while s will be of 
order 1/a, which is large. In the end it will hopefully turn out that the result 
will not depend on what kind of "suitable renormalization scheme" we used, 
however, we are not yet at this point so we should still make sure that we do 
not have large expansion coefficients. 



3 Expansion in 1/Xi 

For an expansion to work, we need to know that we are expanding with respect 
to a small parameter. The idea of perturbation theory is that the coupling 
constant, a, is a small quantity. Equation ||2Jl can be expanded in powers of a. 
We have 



+ ^ log — + - - —J- - — 3^ + 



Poa 'Pi "KPiaJ ' (31 /?2 2/34 2/32 ' /?a 
3/32 + 3/33 f3l + 3/33 + 3/35 ■ 



(12) 



From this series we see that if a is small, s must be large. Indeed, in the case 
of QCD, s is given by 

where fx is the renormalization scale as it occurs in dimensional regularization 
and Aqcd is of the order of the energy range where the coupling is large and 
perturbation theory is not accurate. Therefore, the perturbation series can 
also be written as an expansion in 1/s. For this we use the inverse of the 
expansion El We have 

^J_.J^r .A /^l r2 

Pf , P3 2/33 ^ , 3/32/3i ^ , /33 5/33 . 



where L, = log(-/3i/(/32s)). 

Now we assume that a suitable renormalization scheme has been chosen, and 
hope to obtain results that turn out not to depend on our "suitable renormal- 
ization scheme" and hence might also have been calculated if we had started 
out with an unsuitable renormalization scheme. In a suitable renormalization 
scheme, we expect that expansion coefficients are not large. In particular ri is 
expected to be of order unity. From equation (|11() it follows that Xi = ri+PorQS. 
We conclude that Xi must be a large quantity, because s is. Hence, the expan- 
sion 

is a good expansion. We substitute this equation into equation H14|) obtaining an 
expansion of a in 1 / Xi . This expansion is substituted into equation Q . We then 
obtain an expansion of the physical quantity R into 1/Xi. Furthermore, in this 
expansion, we substitute for r2, r3, r4, . . . the values as given by equation pi|l . 
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It should be noted that although Xi does not depend on the renormalization 
scheme, it is dependent upon the physical quantity under consideration. There 
is, however, nothing wrong with using a different expansion parameter for every 
different physical quantity. The expansion that is obtained by making all these 
substitutions is 



(16) 



AX, , 3^^/^^^^L. + mf-i-2L. + Ui + L3 



where Lx is given by Lx ~ log(— ro/3i/(/3o-'^i))- We see that in this expansion 
all scheme-dependent terms have canceled. These scheme-dependent terms are 
the ones involving r\ and /3i with i > 2. Because of the cancellation of these 
terms, we have obtained a scheme-independent perturbation series. 

For this method to work for any order in perturbation theory, we must 
prove that the cancellation of scheme-dependent terms happens at every order 
and not just up to fourth order. For this purpose, note that our expansion is 
an expansion with respect to the variables l/^i, 1/ Po, Pi, (32, Ps, ■ ■ ■ ,rQ,ri, X2, 
X,, X4, . . . , Lx- If we refer to the order of a term in the series, we mean the 
order in l/Xi. We should prove that actually the variables ri, P2, P3, . . . do not 
occur in this series expansion. Let us assume that such a variable actually does 
occur at some order n in the expansion. 

We introduce v as an alias for one of the offending variables (there might 
be several offending variables) that occurs at order n. This means that dR/dv 
is of order n. We consider what happens if we re-expand the expression for 
i? in a again, and then differentiate with respect to v. To do this we first 
have to expand 1/Xi — l/{PorQS + ri) in 1/s and then use equation H12|l to 
expand this in a again. Xi was defined in such a way that it actually does not 
depend on scheme-dependent quantities such as v. Therefore, we know that 
the entire series of l/Xi in a docs not depend on v. Furthermore, we note 
that during re-expanding and differentiating the order of a term in 1/Xi, 1/s 
or a, whichever applies, never decreases. Therefore, all the invariant terms up 
to order n after differentiation cause terms that are at least of order n + 1. The 
only terms that can give a contribution of order n are the scheme-dependent 
terms. However, the terms obtained by differentiating a, using the chain rule, 
are of higher order, so this does not contribute. The conclusion is that the 
derivative with respect to v up to order n is the same before re-expanding as 
after re-expanding, except that 1/Xi is to be replaced by —ajr^ and Lx is to be 
replaced by log(/3ia//3o). Hence, the physical quantity R up to order n depends 
on the scheme at order n. This is a contradiction with the starting point of 
section 121 Therefore the expansion in must have renormalization scheme 

independent coefficients. 



4 Resumming the Logarithms 

In this section, we will resum all logarithms Lg that occur in equation (|14(l . We 
start out by observations that have been made from computer algebra experi- 
mentation, but in the end we will prove our results to be correct to all orders. 
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By looking at the expansion 114|) . we observe that the leading logarithms, i.e., 
terms of the order i"~^/s" can be summed into the quantity 1/s' defined by 



1 1 1 



s' sl + PiL/{Pls) 
After this resummation we have the expansion 



(17) 



' ' ' ' ' 2L-\L^- \L^) (18) 



V 6 ' 2^ 3 



/?8(s')5^ ' 3/37(s')5 6/37(s')5 ZPlis'f 

Looking at the highest order logarithms in this expansion, we recognize the 
expansion of /3i/(/3oS^) log(l — PiL/{(3qs')). We therefore define a quantity L' 

by 

L' = log fl - . (19) 



It now turns out that after rewriting the expansion for a with respect to 1/s' and 
L', we recover the original expansion where 1/s has been replaced by 1/s' and 
L hy L' . Because of this we can iterate this procedure and obtain a sequence of 
values l/s„ and from the iteration 



/3i 



Ln 

(20) 



L„+i = log 



1 A- -4^T 



where we have rewritten in terms of s„ and i„ instead of in terms of 

Sn+i and Ln- We note that this iteration increases the order of L with respect 
to 1/s. Therefore, in the perturbative regime, the iteration should make L 
converge to zero. This resums all logarithms into a new value for 1/s. 

A curve can be drawn through the sequence of points (l/s„, L„). This curve 
is given by 



111 



S{x) Soo 1 

L{x) = log 



/3i 



1 



(21) 



1 



Pi 



where x parameterizes the curve and different curves (for different initial values 
of 1/s and Lg) are labeled by 1/soo- That this is correct, can be seen by 
checking that the iteration for the pair (1/s, L) is recovered if x is iterated using 
the prescription 

Xn+l = log — . (22) 

We see that also for the a;„ the property holds that Xn+i is of higher order than 
Xm hence in the perturbative regime it should converge to zero. In this case 
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we note that l/s{x) I/sqo, which explains our notation "1/soo" to label the 
different curves. 

We must still prove our assertion that the iteration preserves the expansion 
to all orders in 1/s. It is sufficient to show that the value of a is constant along 
the curves introduced above. The infinitesimal form of the curves is 



5s = ^5x; 

Po 

5Ls = ( -1 - 1 5x 



(23) 



Proving that this is a symmetry of the expansion H14() is also sufficient to see 
that the iteration works. This is what we will do in the next section. 



5 Proof of the Iteration 

Here we will prove that the transformation l|23|) is a symmetry of the expansion 
equation (|14|l . This, at the same time, shows that the iteration of equation H20() 
works to all orders and that the logarithms in the expansion can be made zero. 
We introduce a quantity a that is a power series in 1/s and Lg. This quantity 
has the definition 

da 

da 

[sa(l/s,L, = 0)],_^ 

We will show that this quantity a is actually the same as a. First note that these 
differential equations are consistent because d/{d{l/s)) and d/dLg commute. 
Secondly, if we confine ourselves to the surface Lg — log(— /3i/(/3qs)) we have, 
as is verifiable by using the chain rule for differentiation, da/ds = P{d), so, also 
using the boundary condition, we sec that on this surface a and a are the same. 
We now show that for every order in 1/s, there is a finite number of terms. This 
ensures that also away from the surface Lg = log(— /3i/(/3gs)) these functions 
are the same, because it is impossible to express s and Lg into each others in a 
finite number of terms. The pre-factors of the expansion in 1/s and Ls of a can 
be obtained by setting as — — l//3o in the derivatives 9"+™a/(9L™9(l/s)"). 
At first sight, it may seem that problems could be caused by terms of the form 
a™s" where n > ni. However, the fact that the differential equations have a 
solution that is an expansion in L]p/s" with n > and m > 0, ensures that 
these problematic contributions will cancel if we substitute as — > — l//3o. Still, 
before this substitution is made, there will be terms of the form a^ {1 / s)'^^'' , with 
p > and q > 0. If we consider the quantity d"a/d(l/s)", the maximum value 
of q for which this type of monomial will occur is equal to n — 1. Differentiations 
with respect to Lg increase the order in a, hence, in 9"+'"a/(9L™9(l/s)"), the 
maximum value for q will be n — 1 — m. Hence, if to > rt we will only have 
terms containing a^(l/s)'' with p > q. If we substitute as — — l//?o, taking 
s — > oo, these terms will become zero and thus do not contribute. From this we 
see that the maximum order in Lg that occurs in the coefficient of l/s^ in the 



■ = §m 

Po 

= -l//3o. 



Pis 
PI 



P{d)- 



(24) 
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expansion in s is p. Hence, this cocfBcicnt of l/s^ contains a finite number of 
terms, as we set out to show. We conclude that the quantity a indeed has the 
same expansion with respect to 1/s and Ls as a. 

Furthermore, from equations (|23|l and (|24|l . it foUows that da/dx = 0. There- 
fore the symmetry (I23|l is indeed a symmetry of the expansion of a in 1/s and L^, 
and the iteration in equation (|20|) keeps the value of a constant to all orders 
in 1/s. 



Equation 1211 then gives 



6 Resumming Lx 

The symmetry of equations 123(1 can be turned into a symmetry of the expan- 
sion of a physical quantity in 1/Xi and Lx- This will enable us to perform 
resummation of logarithms in such an expansion. From Xi = ro/?os -I- ^i, it 
follows that 

6Xi - roMSs) = ^Sx. (25) 
Po 

Turning Ls into Lx is done via 

Ls^Lx+\og . (26) 

Xi - ri 

Applying S on both sides gives 

-l-^)sx = SLx-^SX,. (27) 
Pos J Xi 



SLx^-{l + ^)sx. (28) 

Using this symmetry it is possible to turn Lx into zero, thereby ridding ourselves 
of logarithms. The value of Xi that is obtained while turning Lx to zero will 
be called Xi. Integrating equations H25I) and (|28|l . we obtain the equation 

Xi =Xi--— log— (29) 

This can be expressed in the Lambert W^-function. This function is by definition 
the solution to W^(z)e'^(^^ = z. We have 

X, = !:^w(-e^). (30) 
Po V / 

Hence, the conclusion is that we turned the standard perturbation theory into 
an expansion in the quantity 1/Xi. The expansion looks as displayed in equa- 
tion H16|l with Xi replaced by Xi and all terms that have an Lx removed. Since 
this reduces the number of terms considerably, let us display a few more. We 
have 

X 



^ 413! 






1 3ra/32X2 




XI 


2PlXt Xf 6/34xf 


Pix! X! 


17 rlPl 


iriaix^ 


iriPfX, r^,X, 




uplxf 


2[3lXt 


Pixf xf ■ 





(31) 
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highest order 



not using symmetry 



using symmetry 



1/. 

l/s4 



91.5 
282.2 
247.5 
248.6 



91.5 
249.0 
249.3 
249.0 



Table 1: Values for Aqcd in MeV obtained by solving equation H12|l to various 
orders numerically, while using or not using the symmetry in equation H23(l . We 
used a number of quarks equal to five to obtain this result. 



In their Review of Particle Physics |S] the Particle Data Group suggests using 
equation (|14() to define Aqcd- To be fully accurate, their definition is not 
completely the same. Instead of our Lg = log(— /3i/(/3os)) the PDG uses Lg = 
— logs. This amounts to a shift in the parameter s. A way to see this is from 
equation |(5J). Adding a factor — /3i//3o inside the logarithm in this equation, 
turns our expansion into the one of the PDG. This is equivalent to adding a 
constant to the right-hand side of this equation. This constant can then be 
moved to the left-hand side, so we see that s is indeed shifted. The consequence 
is that the Aqcd that we use differs by a multiplicative constant from the one 
of the PDG. We have 



where the beta function coefficients are given in our conventions. In the rest of 
the paper we use our conventions, hence we will be writing Aqcd for AqJ^^^ 
and never mention Aq^§ again. Note however, that the expansion with respect 
to 1/Xi and Lx becomes different if we choose a different prefactor inside the 
logarithm. We will see in the next section that using the symmetry of equations 
(|25|l and H28(l resolves this ambiguity. 

The suggestion of the PDG to use equation H14|l to define Aqcd is not very 
practical. In the first place it would seem to be easier to use equation |(2Jl. 
Secondly, if one is going to use equation H14|l . the symmetry from equation 123|) 
is useful to obtain a series that converges faster by transforming Ls to zero. 
The PDG gives as{Mz) ~ 0.1187. In tableHwe compare the value obtained for 
Aqcd using the symmetry and not using it. The used beta function coefficients 
were calculated in [H] and recently confirmed in . The number of flavours was 
set to five. We indeed see faster convergence. 

Below we will be using Aqcd- The value that we use comes from equa- 
tion H14|l to fourth order where we use our symmetry to get rid of the L^'s. 
If one uses standard perturbation theory, one has to pick a suitable value for 
/I and specify the renormalization scheme. In that case, the consistent way to 
proceed is to only use the beta function up to the loop level to which the rest 
of the calculation is done. The difference with our case is that while in the 
standard approach as{Mz) ~ 0.1187 appears as a fundamental constant, in our 
approach Aqcd would be the fundamental quantity. Determining the value of 
a fundamental constant from another one is better done with as much accuracy 
as possible. This is the reason that we use the beta function up to four loop 
level. 



7 Determining A, 



QCD 




(32) 
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8 How Invariant is Invariant, Really? 



Our method attempts to remedy the arbitrariness in choosing the renormahza- 
tion prescription, so we should now ask the question to what extend this method 
itself is arbitrary. A first possible source of arbitrariness is the choice of variables 
to parameterize the renormalization prescription. We choose the set of param- 
eters ri, /?2, /?3, /34, . . . to parameterize the renormalization prescription. If one 
chooses a different set of variables, fi = ri + Ai, /32 — 02 + '^2, = /33 + A3, 
where the A's are constants, the integration constants Xi from section|5]also be- 
come different. However, if we have Ai = 0, it turns out that in the end the final 
expansion coefficients still have the same value, so this is not an arbitrariness 
of our method. 

This shows that the small invariant quantity that we decide to use for ex- 
pansion is more important than the precise definition of the other invariants. 
The question that might arise is what would happen if we would expand with 
respect to some arbitrary function of Xi instead of with respect to Xi. We 
could, for instance, expand with respect to the sine of Xi. The possibility of 
expressing the coupling constant in one scheme as a power series in the coupling 
constant in another scheme is a possibility that has been mentioned in liter- 
ature, for instance in |10| . Our point of view is that the possibility to use an 
arbitrary power series is not a fundamental arbitrariness of perturbation theory. 
Note that if it were, it would apply to any perturbative method in any branch 
of physics. Small quantities that are used for expansion should be the ones that 
come naturally with the problem under consideration, not the ones that can 
be used to show that any approximation method can be made to give wrong 
answers. In field theory, the situation is that no renormalization scheme is a 
priori better than any other, and this ambiguity can be parameterized by using 
an arbitrary power series in the coupling constant. This is the reason arbitrary 
power series of the coupling constant can be useful to consider. Considering 
the fact that the first scheme invariant arises naturally from the demand that 
a physical quantity should not depend on the scale, it does not make sense to 
consider arbitrary functions of this, perhaps with the exception of a translation 
in de definition of Xi. I.e., a non-zero value of Ai. 

If Ai 7^ 0, we obtain a different invariant variable on the one loop level, 
namely X( = Xi -|- Ai. If we were to expand in l/X[, we would indeed obtain 
a different expansion. Here, it appears, we have finally found arbitrariness. 
However, we actually already mentioned this. An example of this arbitrariness is 
given in the paragraph that contains equation (|32(l . The prefactor that we choose 
inside the logarithm Ls is equivalent to choosing a value for Ai. Resumming 
the logarithms yields the same result for Xi and, hence, also the same result 
for the physical quantity R. 



9 Hadronic-i? 

In this section we consider massless QCD-corrections to hadronic i?had- This 
quantity is by definition given by 

cr(e+ + e- hadrons) 

^had - — T-j— Z 1—, TT' ['^'^ 

a{e+ + e ^ 11+ + n ) 
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Figure 1: hadronic i? as a function of the center of mass energy for Nf = 5, as 
found by various methods, normalized by dividing by the PDG-one-loop result. 
Note that our method gives rather different results from the PDG-method and 
the PMS-method. From top to bottom we see graphs for PMS-two-loop, PDG- 
two-loop, [gap], our- two- loop, PDG-three-loop, PMS-three-loop, [gap], and our- 
three-loop. 



where electroweak corrections are neglected. The value of this quantity can 
be obtained up to third order from the review paper The beta function 
coefficients can be obtained from the same paper or from [H]. a is taken to be 
otsji^. For the first renormalization scheme invariant Xi, we have 

Xi = ri + /Joros 

.-_„«3,-»,(ii-|c,3,).(-H,i..),.(^), 

(34) 

where we have written scm for the squared center of mass energy, to avoid 
confusion with the quantity log(/i^/AQ(-.j3), that we also call s. The addition of 
(3orQS has resulted in the replacement fi Aqcd, in ri. Note that the value of 
Aqcd is scheme dependent. For instance, if we would use the MS-scheme (we 
are actually using MS), the two values of Aqcd would differ by a multiplicative 
constant. The scheme invariant is independent of this choice. A practical way 
to proceed is to start from the expansion of -Rhad up to and including terms of 
order a^, substitute, equation (|14|l into it and expand up to terms of order 
In this we substitute equation (|15l) and expand again with respect to 1 /Xi up 
to terms of order l/Xf. Using computer algebra this is a more or less trivial 
thing to do. This cancels all dependence on the scheme, as we have argued. 
In particular, we observe that results no longer depend on the renormalization 
mass fi. 

We plot a normalized variant of hadronic R as calculated by two other 
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Figure 2: hadronic i? as a function of the center of mass energy for Nf = 5, as 
found by various methods, normalized by dividing by the PDG-one-loop result. 
As opposed to the previous graph, our method now gives very similar results to 
the PDG-method and the PMS-method. From top to bottom we see graphs for 
improved-our- two- loop, PMS-two-loop, PDG-two-loop, [gap], PDG-three-loop, 
PMS-three-loop, improved-our-three-loop. 



methods that have been proposed to handle the scheme dependence. The first 
method is the one used by the Particle Data Group. The PDG uses MS and sets 
A* — \/scM- The second method that we consider is the PMS-criterion. Infor- 
mation on this can be obtained from IQ . We also plot the result obtained by 
our method of re-expanding in 1/Xi. We obtain figure ^ In this figure we did 
not yet use the symmetry of equations (|25|l and H28() to remove the logarithms. 
If we use it, we get figure |2 The normalization that we mentioned is done by 
dividing by the one-loop result that can be found using the PDG-method. This 
prevents the graphs from being very close to each others. We conclude that our 
method gives results basically equal to the two other methods provided that 
it is improved by the use of the symmetry. "Basically equal" means that the 
uncertainty that results from ignoring the next order is much larger than the 
uncertainty that comes from the renormalization scheme dependence. 

10 Comparison with Other Solutions 

Other solutions to the problem of renormalization scheme dependence have been 
proposed. The one that perhaps is most like ours, is by C.J. Maxwell In 
fact, in the case where /3i — the one- loop result of Maxwell is identical to 
ours. He sums some of the higher order terms along with the ones that come 
from orders where the full result is known. If an n-loop calculation has been 
performed, Xn is known and all terms that contain Xi with i < n are to be 
summed. This has the problem that the definition of the Xi depends on what 
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parameters are chosen to parameterize the scheme. Maxwell notes this himself 
in [H]. As we have seen, the invariants are constants of integration, hence they 
depend on the boundary chosen. Therefore it does not make much sense to 
resum these, because it will inevitably lead to arbitrariness. Our approach of 
expanding in the first invariant avoids this problem. 

The PMS criterion, introduced in |[10) has the big advantage that, apart from 
our own method, it is the only one that is really completely scheme and conven- 
tion independent. It has, however, the disadvantage that it is difficult to apply. 
Here, the optimum value, in some sense of optimal, of scheme dependent pa- 
rameters is determined. This involves solving trancedental equations containing 
integrals, and hence can generally only be done numerically. If one is interested 
in expressing a physical quantity in, say, the numbers of flavours Nj, this can 
only be given as a set of equations. By contrast, we have a series expansion 
in the parameter 1/Xi. Only one equation needs to be solved to obtain Xi 
from Xi. Such a result can easily be expressed in, for instance, Nf. 

The method of calculating hadronic R used by the Particle Data Group 
consists of setting /i to a "good" value. This has the disadvantage that one has 
to pick this "good" value alongside with a "good" renormalization scheme. For 
hadronic R the PDG makes /z equal to the center of mass energy and uses MS. 

Another idea was put forward in jHj. Here we obtain a differential equation 
for the physical quantity. We have 

f{R). (35) 



dscM 

As before, sqm is the energy squared in the center of mass, not the scale. It 
turns out that the right-hand side is scheme independent, which should not come 
as a surprise to the reader of this paper. The method has some disadvantages 
compared to ours. Firstly, this differential equation still needs to be solved. 
This could be done by giving an initial value at some reference energy, and then 
expanding in log(scM/sref)- This, of course, goes wrong if the energy of an 
experiment starts to differ significantly from the reference energy. Furthermore, 
the mere mention of a "reference energy" indicates that we are reintroducing 
arbitrariness. So, presumably, we are not to solve this differential equation by 
expanding with respect to this quantity. What are we to do then? Solving the 
differential equation numerically, perhaps? The reader will not find it difficult 
to think of disadvantages of this. Furthermore, we would want to relate different 
physical quantities to each other. In this method we would give a series expan- 
sion that expresses one into the other. However, which of all possible physical 
quantities is going to appear in a listing of fundamental quantities? And at 
what energy is this quantity going to be listed? This method has no preference 
for a particular quantity or energy. 

Yet another idea can be found in |3]. Here an equation that looks a lot like 
equation |(5J is given. In our notation it would be given by 

^ , ro^, roPo + PiR f 1 1 

R /3o f3iR Jo \Y{R') YoiR'Y +Y,{R'f 

(36) 

where Y{R) is a power series that starts with the term Yol-R)^- The Yi are 
invariants related to the Xi but their definition is not entirely the same. If 
we go to a scheme where all with i > are zero, the just given equation 
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and equation |(2J) are the same. Because only renormalization scheme invariants 
occur in eauation l36l it must hold in any scheme. This could be an answer to the 
question that bothered us in the previous paragraph, because this expression 
is a solution to the differential equation (|35() . The relation to our method is 
that we propose to simplify this expression by turning it into an expansion. 
Numerical accuracy is then achieved by resummation of logarithms. 

The method proposed by @] ultimately boils down to the same equation . 
The philosophy is different though. The idea is that the coupling constant 
should be interpreted as an "effective charge" that no longer receives higher 
order corrections. 

11 Conclusions 

We can get rid of the unphysical dependence of physical quantities on the 
renormalization scheme by expanding in 1/Xi where Xi, is the renormaliza- 
tion scheme invariant that occurs on the one loop level. Numerical accuracy is 
achieved by resumming all logarithms that contain Xi into the quantity 1/Xi. 

For future research it would be interesting to look at the possibility of gen- 
eralizing our approach to the case of a theory with masses and/or multiple cou- 
pling constants. Another thing that could be done is to try a similar approach 
to factorization scheme dependence as it arises when one studies deep-inelastic 
scattering. In the context of the method introduced in [S], both renormalzation 
and factorization scale dependence are discussed in [7]. 
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